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ABSTRACT 



The ion plays a key role in the chemistry of dense interstellar gas clouds where stars and planets are forming. The low temperatures 
and high extinctions of such clouds make direct observations of impossible, but lead to large abundances of H 2 D + and D 2 H + , which 
are very useful probes of the early stages of star and planet formation. The ground-state rotational orfho-D 2 H + lij-Oo,o transition at 
1476.6 GHz in the prestellar core 16293E has been searched for with the Herschel** HIFI instrument, within the CHESS (Chemical 
HErschel Surveys of Star forming regions) Key Program. The line has not been detected at the 21 mK km s -1 level (3cr integrated 
line intensity). We used the orfho-H 2 D + li,o-li,i transition and para-D2H + li,o-lo,i transition detected in this source to determine an 
upper limit on the ortho-to-para D 2 H + ratio as well as the para-D 2 H + /orfho-H 2 D + ratio from a non-LTE analysis. The comparison 
between our chemical modeling and the observations suggests that the CO depletion must be high (larger than 100), with a density 
between 5 x 10 5 and 10 6 cirr 3 . Also the upper limit on the orfho-D 2 H + line is consistent with a low gas temperature (~ 11 K) with 
a ortho-to-para ratio of 6 to 9, i.e. 2 to 3 times higher than the value estimated from the chemical modeling, making it impossible to 
detect this high frequency transition with the present state of the art receivers. 

Key words, astrochemistry - ISM: individual (16293E) - ISM: abundances - Line: identification - Radiative transfer 



1. Introduction 

In the recent years, the chemistry of dark clouds and star forming 
regions has constantly been revise d with the discovery of multi 



ly deuterated molecules: D2CO dTurneij[ l990: Ceccarelli et al. 
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1 998 k ND^H fRoueff et al.ll200(ih. D,S (IVastel et al.l2003l) . ND 
dLis et alJ 120021: Ivan der Tak et alj 120021). CHD 2 OH, C D3OH 



dParise et al.l I2002L 120041) and DoH + (IVastel et alJl2004l) . Two 
main pathways can be invoked for understanding the observed 
large deuterium fra c tionation. The first is ba sed on grain chem- 
istry (lTielenslll983h . lVastel et al.l (l2003h and lParise etail {2004) 
showed that grain chemistry models require a very high atomic 
D/H ratio accreting on the grains in order to explain their high 
deuterium fractionation ratio that could not be reproduced by 
gas-phase modeling at that time. New models can now reproduce 
the observed abundance of the multi-deuterated isotopologues 
of formaldehyde an d methanol, both f ormed in the last stage of 
the prestellar phase (Taquet et al. 2012). The second pathway for 
forming deuterated molecules is based on gas-phase chemistry 



* The chemical network is available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A 

** Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



and results from the ion-molecule deuterium exchange reactions 
taking place at low temperatures. In this scenario, deuterated 
molecules are produced through successive reactions starting 
with H2D + , dominant at temperatures lower t han 20 K, CH2D + 
or C? HD + , dominant at higher temperatures dRoberts & M illar 
2000; iGerlich et afll2002l) . H 2 D + has proven to be a very good 
probe of the dense cold cores, where CO disa ppears from the gas 
phase and is depleted on t o the dust grains (ICaselli et al . 2003, 
120081) . iPhillips & Vastell d2003l) pointed out that the deutera- 
tion should be extended beyond H2D + , to D2H + and and 
suggested that the detection of the D2FP ion might be possi- 
ble. Calculations including all possible deuterated isotopomers 
of have confirmed that, in dense CO depleted regions, the 
abundance of D2I-P will be similar to that of HyD + , and the 
D/H ratio will be largely enhanced dRoberts et alj|2003l [2004; 
IWalmslev et alJl2004t lCeccarelli & Dominikll2005l) . 
D2I-P in its para form w as detected f or the first time in the 
prestellar core 16293E by IVastel et al.1 (120041) . showing its im- 
portance, as well as that of D3 , in determining the total deu- 
terium abundance i n the gas phase. Th is dense core, revealed by 
amonia emission (Mizu no et al.ll 19901). is sh eltered in the dense 
cloud L1689N (see IWootten & Loren1ll987h . in the Ophiuchi 



region (distance 120 pc: Loinar d et al 



2008) and has been re- 



vealed by millimeter lines and continuum emission. This cloud, 
also harboring a young binary protostellar object (IRAS 16293- 
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2422 A and B) and bipolar outflows (named Rho Oph East 
by iFukui et aD (119861) '). has been extensively observed, reveal- 
ing an extreme molecular deuteration in particular t oward s the 
cold dense core, named 16293E by lLoinard et all d200lh and 
ICastets et alj (|200jJ). Note that this core has been called p Oph 
E (as in Rho Oph Eas t referring to the o utfl ow emanating from 
IRAS 16293-24221 bv lSaito et al.1 (l2000h and lGerin et all (1200 lh . 
However, we decided to use the name 16293E since p Oph E 
already refers to another c ondensation withi n the L1688 cloud 
in the Ophiuchi complex dLoren et al.l 1 1990h and, therefore, is 
confusing. We present here recent Herschel/ HIFI observations 
obtained within the Herschel guaranteed time Key Program 
CHESS, and modeling of both ground state transitions of para- 
and ortho-D2H+, as well as ortho-H2D + , us ing the recent col li- 
sional coefficients with ortho- and para-H 2 dHugo et alj|2009h . 

Figure Q] shows the 1.3 mm continuum map as well as the 
DCO + contour map (blue) dLis et al.ll2002l) . Studies of the gas 
kinematics using CO, HCN, H 13 CO + , HCO + and DCO + trac- 
ers lead to the conclusion that the deuterium peak is a part of 
the a mbient cloud t hat is pushed and compressed by the out- 
flow dLis et al. 1 12001 . Indeed this shock could have released 
deuterated species, that were condensed on the dust grains, into 
the gas phase with a subsequent cooling of the gas to lower 
temperatures. This effect, combined with the low temperature 
gas-phase chemistry in the high-density shock-compressed gas 
leads to a high mo lecular deuteration as obs erved in: PNC (9% : 
Hirota et al.ll200lh. D 7 CO (4 0% + 20% lLoinard et all l200lt 
Ceccarelliet al.ll2002|). N?D + 



200 lh. NH^D (19% 



_ - dGerin et al.l 

Mizuno et al l Il990t lLoinard et alj 120011: IRoueff et alj 120051). 
ND 2 H (4%: iMizuno et alJll990t lLoinard et alJl2001t iLis et alJ 
2006) , DCO + an d DCN (~ 10%: ILis et alJ 120021). HPS 
dVastel et al.ll2003l) . D 2 H + dVastel et al Il2004h . HDO dStark etal] 
|2004j), ND 3 (0.1%: IRoueff et alj|2005h . Because no far-infrared 
or submillimeter point source has been found, the source can be 
classified as a prestellar core, the early phase of the formation of 
a protostellar object, before gravitational contraction occurs. 



2. Observations and data reduction 

The para-D 2 H + ground state transition was observed to- 
wards the prestellar core 16293E with the HIFI instru- 
ment (de Gra auw et al. 20101) o n board the Herschel Space 
Observatory dPilbratt et al.l 201 Oh . as part of the Hersche l guar- 
anteed time Key Program CHESS dCeccarelli et alJl2010h . 

A pointing at a frequency centered on the D 2 H + line (~ 
1476.6 GHz) with the band 6a HEB receiver was performed 
on February 16, 2011 using the pointed Double Beam Switch 
(DBS) mode with optimization of the continuum. In this 
mode, both the HIFI Wide Band Spectrometer (WBS) — 
providing a spectral resolution of 1.1 MHz (~ 0.2 km s~') 
over an instantaneous bandwidth of 2.4 GHz — and the HIFI 
High Resolution Spectrometer (HRS) — providing a spectral 
resolution of 125 kHz over an instantaneous bandwidth of 0.12 
GHz — were used. The DBS reference positions were situated 
approximately 3' East and West of the source. The HIFI beam 
size at the observed frequency is about 14", and the main beam 
and forward efficiencie s are about 0.72 and 0.96, respectively 
(Roelfsem a et al.ll2012h . The data have been processed using the 
standard HIFI pipe line up to level 2 with the ESA-supported 
package HIPE 8.0 dOtdl2010D . The on-source integration time 
for this observation was 18157 seconds. FITS files from level 2 
data were created and translated into CLASS/GILDAS format 
for subsequent data reduction and analysis. For each scan, 
a low order polynomial baseline was fitted outside the line 




Fig. 1. Continuum map at 1.3 mm convolved to 15 with black 
contours at 4%, 6%, 8%, 15%, 30%, 60%, and 90% of the 
peak located at the IRAS 16293-2422 position (6.7 Jy in a 1500 
beam). Dashed blue lines show the distribution of the DCO + 3-2 
integrated intensity between 2.8 and 5 km s _1 , with contour lev- 
els contour levels are 50%, 60%, 70%, 80% and 90% of the peak 
(integrated over velocities between 2.8 and 5 km s~'). The cross 
represent the position towards which H 2 D + and D 2 H + have been 
observed. 

16293E 
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Fig. 2. para-D 2 H + (1 l,o _ lo,i) an d ortho-D 2 H + (1 i,o-lo,i) transi- 
tions observed with the CSO and Herschel/HIFI. The Wide Band 
Spectrometer is represented in red, and the High Resolution 
Spectrometer in black. 



window. The antenna temperatures were finally converted to the 
T m b scale, using the theoretical values of the main beam and 
forward efficiencies given above. We present in Figure [2] this 
observation as well as the para ground state transition observed 
at the Caltech Submillimeter Observatory (CSO) (Vasteletal. 
2004). Both observations were pointed towards the DCO + peak 
emission, at coordinates crp ooo = I6 h 32 m 28 s .6 2, 5 2 ooo = _ 24° 
29' 2.7" (see for example Roue ff et al.l 120051 Figure 9). Note 
that the coordinates in IVastel et aU (|2004) were not correctly 
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quoted in the text and should be replaced by the here quoted 
coordinates. 

The two nuclear-spin species of H2D + and D2H + (ortho and 
para) are considered separately. The parameters of the ortho- 
and para-D2H + and ortho-H2D + ground transitions are obtained 
using the CASSIS0 software, which takes into account the or- 
tho and para forms separately, with an independent computa- 
tion of the partition function (cf. Formalism for the CASSIS 
software, http://cassis.irap.omp.eu/), and are reported in Table 
1. The Einstein coefficient for the ortho-D2H + transition (~ 3 
10~ 3 s _1 ) is much larger than for the para-D2H + (~ 4.6 10~ 4 
s -1 ) and ortho-H2D + (~ 1.1 1CT 4 s -1 ) transition, since it varies 
as a function o f v 3 . Note that w e ado pt here the most recent 
measurement of lAmano & Hiraol d2005l) of the para-D2H + line 
frequency (691.660483 GHz), 372.421385 GHz for the ortho- 
H 2 D + transition and 1476.605708 GHz for the ortho-D 2 H + tran- 
sition dAsvanv et al.ll2.008b . From a simple Gaussian fitting func- 
tion, the resulting V/ sr are (3.78 + 0.02) km s _1 for the para- 
D 2 H + transition and (3.59 ± 0.02) km s" 1 for the ortho-H 2 D + 
transition (statistical uncertainty from the fit). The minimum dif- 
ference between the lines is 0.15 km s _1 , compatible with the ~ 
0.1 km s _1 resolution of the acousto-optic spectrometer at 692 
GHz, and 0.04 km s at 372 GHz. A difference in the observed 
Vi sr i s found (exploitation of the data mentioned in iLis et alJ 
2002, but not published) in the DCO + 5-4 (3.65 ± 0.01 km s" 1 ) 
and 3-2 (3.55 + 0.01 km s _1 ) transitions. The quoted uncer- 
tainties correspond to the Gaussian fit and the line parameters 
are listed in Table 1. A larger difference in the line center ve- 
locity is noticeable comparing the ortho-ND2H li,i-Orj,o (V/ Jr 
- 3.64 km s" 1 ) and 1 1,1— Oo,o (V; sr = 3.60 km s -1 ) transitions 
dLis etal J 120061) . as well as the ortho-NH 2 D li.i— lo.i (V/ Jr = 
3.65 k m s" 1 ) and para-NH 2 D 1 1 . i — 1 o, l (V /jr = 3.67 km s -1 ) tran- 
sitions (lLoinard etal] 1200 ll). with the ND 3 l -0 (Vi sr = 3.35 
km s -1 ) transition dRoueff et alj|2005l) . Such a difference could 
be explained by dynamical motions in this source that cannot be 
reproduced with a simple Gaussian fitting function. The inter- 
action between the outflow of IRAS 16293-2422 and the cold 
core 16293E is likely to produce a velocity shift compared to 
the L1689N parental cloud, with a difference between the non 
deuterated species and their increasingly deuterated forms. The 
difference in the beam sizes of the numerous observations may 
be responsible for such a variation as well. The single-dish ob- 
servations do not allow to conclude. 



The 3cr upper limit on the integrated line intensity, using the 
Wide Band Spectrometer, has been derived following the rela- 
tion: 



3cr (K km s" 1 ) = 3xr«wx V2 x dv x FWHM 



(1) 



with rms (root mean square) in K, dv, the channel width, in 
kms~' and FWHM (Full Width at Half Maximum) in kms -1 . 
We assumed FWHM = 0.29 kms 1 , which is the para-D2H + 
emission line width. 



3. Discussion 

3.1. Line widths 

The line widths are listed in Table Q] The expected thermal line 
width varies with the kinetic temperature (T) as: 



Av T (km s~ l ) = 2 V2Zn2 x J—, 
V m 



(2) 



where k is the Boltzman constant and m is the molecular weight. 
The ortho-H2D+ line width corresponds to a 1 1 .3 K gas, compa- 
rable with the 9.2 K temperature found for para-D2H + . However, 
the uncertainty in the fit of the line, using the Levenberg- 
Marquardt method, should be taken into account leading a tem- 
perature of 11.3 + 1.8 K for H 2 D+ and 9.2 + 3.2 K for D 2 H+. 
The resulting range for the kinetic temperature is therefore [9.5 
- 12.4] K taking into account both ions. 

The observed (o-p) ND2 H line w i dth (about 0.45 km s _I for 
the li.i— 0q.q trans itions; ILis et all I20061) . ND 3 (0.44 km s" 1 ; 
iRoueff et al.ll2005h as well as the DCO + line width (about 0.5 
km s _1 fo r the 5-4 transi tion, less optically thick than the 3-2 
transition: ILis et al.ll2002l) and the N?D + line width (0.32 km s -1 
for the 3-2 transition: lGerin et al.l200ll) are systematically larger 
than the thermal line widths of 0.16, 0.12 and 0.16 km s _1 re- 
spectively for a 10 K kinetic temperature. Considering the high 
critical densities for the (o-p) ND2H and ND3 transitions (larger 
than 10 6 cm -3 ), it seems that systematic motions occur even 
in the dense part of the cloud. Also, at large enough densi- 
ties/depletions, one expect to have the light ions left as the main 
tracers of the gas. This further confirms that H2D + and D2H + re- 
main the only tracers of the cold, dense and CO/N2 depleted cen- 
tral region where other molecular tracers are largely condensed 
onto the dust grains, rather than the lower density envelope when 
turbulence takes over. 



100-r 
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691.7 GHz 



l io 



Fig. 3. Diagram of the lowest energy levels (in Kelvins) of the 
D2H + molecule. 



1 CASSIS (http://cassis.irap.omp.eu/) has been developed by IRAP- 
UPS/CNRS. 



3.2. H 2 D + and D 2 H + non-LTE modeling 

The H2 density is a critical parameter for the interpretation of 
our deuterated ions observations. CO observations of the 3-2, 
4-3, 6-5 transitions give a lo wer limit for th e molecular hydro- 
gen dens i ty of 5 x 10 3 cirT 3 dLis et al.ll2002l) . In their analysis, 
ILis et al.l d2002l) used a kinetic temperatu re of 12 K (based on 
NH3 observations of iMenten et al.|[l987l) and a density of 5 x 
10 5 cm 3 , consistent with their observed DCO + (5-4)/(3-2) and 
N?D + (4-3)/(3-2) line r atios. From 450 and 850 pun continuum 
maps Star k et al.l (12004) inferred an isothermal dust temperature 
T = 16 K and peak density of 1.6 x 10 6 cirT 3 . Note however 
that their dust peak emission does not correspond to the peak of 
deuterated molecules (about 15" away) that we are studying in 
the present paper. This could be due to the release of the deuter- 
ated species that were condensed on the dust grains by the in- 
teraction between the outflow from IRAS 16293-2422 and the 
cold core 16293E. This compression is then likely to cool the 
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Table 1. Derived para meters of the orth o-H2D + transitions as well as the ortho and para D2H + fundamental lines. D CO + 3-2 and 
5-4 transitions used in iLis et all (120021) are also quoted using the estimated frequencies from Case lli & Dord (l2005h . Note that a 
discussion on the Vlsr can be found at the end of Section 2. 



Species 


Transition 


Frequency 


Telescope 


beamsize 


J T„„, dV 


rms 


binsize 


AV 








GHz 




(") 


(mK km s~') 


mK 


(kms- 1 ) 


(km s" 1 ) 


(km s' 1 ) 


ortho-H 2 D + 


li,o _ li,i 


372.421385 


cso 


20 


720 


80 


0.080 


0.36 ± 0.04 


3.59 ± 0.02 


para-D2H + 


1 i,o— lo,i 


691.660483 


cso 


11 


183 


170 


0.080 


0.29 ± 0.05 


3.78 ± 0.02 


ortho-D 2 H + 


1 i,i— Oo,o 


1476.605708 


HIFI/HRS 


14 


<30 


120 


0.012 






ortho-D 2 H + 


1 i.i-Oo,o 


1476.605708 


HIFI/WBS 


14 


< 21 


21 


0.203 






DCO + 


3-2 


216.1125822 


CSO 


35 


3235 


94 


0.134 


0.67 ± 0.03 


3.55 ±0.01 


DCO+ 


5-4 


360.1697783 


CSO 


20 


528 


109 


0.081 


0.54 ± 0.03 


3.65 ± 0.01 



gas to lower temperatures. These observations possibly trace the 
parental cloud of the prestellar core 16293E. We will therefore 
consider in the following study densities between 10 5 and 10 6 
crrT 3 and kinetic temperatures between 9 and 16 K (see section 
ED- 

From the recent computations of the collision rates (Hug o et all 
2009), the critical densities for the ortho-H2D + , ortho-D2H + 
and para-D2H + transitions are ~ 1.3 x 10 5 cm~ 3 , ~ 7 x 10 6 
cm -3 and ~ 5. 6 x 10 5 cm' 3 respectively, us ing the Einstein co- 
efficients from lRamanlal & Tennyson] (12004 . 1.2 x 10" 4 , 3.3 x 
10~ 3 and 5.1 x 10~ 4 s respectively. From the large difference 
between the critical densities, the ortho-D2H + transition will 
trace denser regions than the para transition. Since the typical 
density in the core seems to be less than these critical densi- 
ties, a Local T hermodynamic Eq uilibrium approximation is not 
applicable. In lVastel et alj d2004t) . only LTE modeling could be 
performed because the collision rates were not available at that 
time. We produced collisional files, assuming a simple 2 level 
system, using the inelastic state-to-state rate coefficients for the 
ortho and para ground transitions of H2D + and D2FL in collision 
with para and ortho H2, as a function of tem perature. We used 
the no n-LTE radiative transfer code RADEX (Ivan der Tak et alj 
120071) in the large velocity gradient (LVG) approximation with 
the collisional files that were created for these species. We could 
consider in the following non-LTE modeling two values for the 
H2 ortho to para ratio of 3 (highest value at thermodynamic equi- 
librium) and (only collisions with para H2). However, the sim- 
ilarities of the collision coefficients with ortho and para H2 for 
(o,p)H2D + or (o,p)D2H + will lead to a very small difference in 
the derived column densities. 

The variation of the para-D2H + to ortho-H2D + ratio as a func- 
tion of the gas temperature, for densities of 1 x 10 5 (green), 5 
x 10 5 (blue) and 1 x 10 6 (red) cm" 3 , are presented in Figure|4] 
(thick lines) and the upper limit on the ortho-to-para D2H + ratio 
is presented in Figure[5](thick lines). Note that a selection on the 
optical depth has been performed (r < 100), since the excitation 
temperature found using RADEX may then not be representa- 
tive of the emitting region. The optical depth for the ortho-D2H + 
transition is lower than 100 for kinetic temperatures larger than 
~ 10.5 K. 

3.3. Full chemical modeling 

The chemical network used here has been collected from various 
sources. First, a reduced chemical network based on the Nahoon 
code was obtained from the KIDA0 database (Wakelam et al. 
[2012b . The network contains species with three or less elements 
and includes hydrogen, helium, carbon, nitrogen and oxygen. 



2 http://kida.obs. u-bordeaux 1 .fr/models 




n h- . : . . . , : . , , . . 

9 10 11 12 13 14 15 16 



T (K) 

Fig. 4. Variation of the para D2H + to ortho H2D + ratio as a func- 
tion of the gas temperature, for densities of 1 x 10 s (green), 5 x 
10 5 (blue) and 1 x 10 6 (red) cm 4 . The thick lines show the ratio 
derived from observations, using non-LTE modeling, whereas 
the thin lines show the chemical model results. Three depletion 
factors at steady-state are presented: complete depletion (CD: 
solid lines), depletion factor =100 (dot-dashed lines) and de- 
pletion factor = 10 (dashed lines). The thin dashed vertical lines 
present the [9.5 - 12.4] K range. 

The reduced Nahoon code has then been modified to include 
deuterium and deuterated species, with H2, HD and D2 form- 
ing on grain surfaces. Finally, the spin states of all isotopo- 
logues of H2 and H^ were incl uded. The ra t e coefficients of the 
new reactions were taken fr om lSipila et alJ (120101) . iH ugo et al. 
(l2009h . lFlower et alJ (12004 and lWalmslev et all (12004 . choos- 
ing the most recent values in case of multiple choice. The re- 
action balancing routine within the Nahoon program (the one 
which checks the reaction list, to make sure that reactants and 
products have the same number of elements and charges) had 
to be modified after the inclusion of the different spin states, to 
avoid to balance spin states between reactants and products (as 
spin states are not conserved). For this purpose, different spin 
states were labelled adding an extra column in the element and 
species definition file. In this column we assigned -1 to the lower 
spin state and + 1 to the higher spin state of all isotopologues of 
H2 and H3 (para-D+ has been assigned a value of +2 as it is the 
highest of three spin states). All other species have a zero in the 
corresponding column, indicating that no spin state is consid- 
ered. The H2 self shielding data has not been modified but new 
parameters for the self shielding of HD and D2 have been intro- 
duced and have been each set to initially be equal to that of H2. 
Coulomb focusing was taken into account for reactions involv- 



4 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 




Fig. 5. Variation of the ortho-to-para D2H + ratio as a function of 
the gas temperature, for densities of 1 x 10 5 (green), 5 x 10 5 
(blue) and 1 x 10 6 (red) cm 4 . The thick lines show the ratio 
derived from observations, using non-LTE modeling, and can 
be considered as upper limits, whereas the thin lines show the 
chemical model results. Three depletion factors at steady-state 
are presented: complete depletion (CD: solid lines), depletion 
factor =100 (dot-dashed lines) and depletion factor =10 (dashed 
lines). The thin dashed vertical lines present the [9.5 - 12.4] K 
range. 




i i 

green = 1e5 cm-3 
blue = 5e5 cm-3 
red = 1e6 cm-3 



14 15 16 



Fig. 6. Prediction of the ortho-to-para H2D + ratio as a function 
of the gas temperature, for densities of 1 x 10 s (green), 5 x 10 s 
(blue) and 1 x 10 6 (red) cm 4 . Three depletion factors at steady- 
state are presented: complete depletion (CD: solid lines), deple- 
tion factor =100 (dot-dashed lines) and depletion factor =10 
(dashed lines). The thin dashed vertical lines present the [9.5 - 
12.4] K range. 



ing n egatively charged ions on neutral grains (Drai ne & Sutinl 
1 1 987b . Only neutral and negatively charged grains are present 
in the network. The final chemical network includes over 
3,500 reactions involving nearly 130 different species. The 
full chemical network used in this work is available at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A4-A/ Comparison 
with previous work and a parameter space exploration will be 
presented by Kong et al. (in prep.). 



A range of H 2 volume densities (1 x 10 5 , 5 x 10 5 , 1 x 
10 6 cirT 3 ), kinetic temperatures (9, 11, 13, 15, and 16 K), and 
depletion factors (= undepleted abundance / depleted abundance 
= 10, 100, infinity) of all the elements heavier than helium (C, 
N and O) have been explored and followed in time until the sys- 
tem reaches equilibrium. Since the CO molecule is the main de- 
stroyer of the ion (and its deuterated counterparts), we sim- 
ulate three depletion values at steady-state: complete depletion, 
and elemental depletion factor of 10 and 100 corresponding to 
CO depletion f actors of about 1 4 and 140, respectively, at equi- 
librium. From lEis et al. 1 (12001 . C 18 O observations at the deu- 
terium peak of 16293E lead to a CO depletion factor (average 
along the line of sight) of ~10, representing a lower limit for the 
central core traced by D2H + . Considering the overall parameters 
in the chemical modeling we adopt the conservative values of 3 
x 10~ 17 s _1 for the cosmic ionization rate, 0.01 for the dust to 
gas ratio, and 0.1 yum for the grain radius. 
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Fig. 7. Variation of the ortho-to-para H2 ratio as a function of 
the kinetic temperature, for densities of 1 x 10 5 (green), 5 x 
10 s (blue) and 1 x 10 6 (red) cm 4 . Three depletion factors at 
steady-state are presented: complete depletion (CD: solid lines), 
depletion factor =100 (dot-dashed lines) and depletion factor = 
10 (dashed lines). The thin dashed vertical lines present the [9.5 
- 12.4] K range. 

Figure 4-6 show the predictions (thin lines) of the model 
at equilibrium (reached in 2 10 5 years for the complete deple- 
tion case, 5 10 5 years for a depletion factor of 100 and 2 10 6 
years for a depletion factor of 10) for temperatures between 9 
and 16 K for three densities. The para-D2H + /ortho-H2D + ratio 
increases with the depletion factor: indeed the disappearance of 
the CO molecule from the gas phase leads to the reactions of HD 
with with an increasing pr oduction of its deuterated counter- 
parts (Phillies & Vastel 2003). Also the increase of the ortho-to- 
para D2I-F ratio with the depletion factor can be explained by 
the fact that the reactions with HD, ortho-D2 and para-D2 con- 
verting para-D2H + in to ortho-D7H + ar e dominant compared to 
the reverse reactions dHugo et al.ll2.009l) . From equations 4 and 
5 we can say that a low ortho-H2 value leads to a very high 
D3 abundance for densities larger than 10 5 cm 4 . Indeed the 
D3 ion is likely to be th e dominant ion in su ch high-density re- 
gions: e.g. Figure 3 from lFlower e t al. (2006) and Figure 2 from 
Sinil a et alJ (l2010h . The ortho-to-para D2FF ratio consequently 
increases as the depletion factor increases, as seen in Figure [5] 
The ortho-to-para H2 ratio should be taken into account in the 
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overall chemistry, since it is critical for the degree of deuteration 
of H3 . Increasing the depletion factor leads to a local decrease 
of the CO abundance. Therefore H + and Hi will mainly react 
with ortho-H2 (at the same temperature the backward reaction 
is negligible), and not CO. This will lead to a high ortho-to-para 
H2 conversion, consequently a lower ortho-to-para H2 ratio as 
seen in Figure [7] Our chemical modeling also predicts a varia- 
tion of the steady-state ortho-to-para H2 ratio for the complete 
depletion case, a depletion factor of 100 and a depletion factor 
of 10, between 9 and 16 K (see Figure|7). Because of the large 
internal energy (~ 170 K) of the lowest ortho-H2 level (J = 1) 
compared to the temperature range explored in this source, the 
ortho-H2 form is a limiting factor for deuteration. It overcomes 
the energy barrier, leading to ex othermic (i.e fast a nd tempera- 
ture independent) reactions (e.g. lGeriich et al. 20021): 

o - H 2 D + + o - H 2 — > (o, p) - H3 + HD (3) 

p - D 2 H + + o - H 2 — > p - H 2 D + + HD (4) 

The excited nuclear spin state of the D3 ion is removed prefer- 
entially by ortho-H2 in the following endothermic (by only 18 
K) reaction: 

m - D+ + o - H 2 — > o - D 2 H + + HD (5) 

Following iHugo et all (|2009) we have assigned meta-D^ with 
the modification having the lowest ground state energy, corre- 
sponding to the Al representation of the symmetry group S3. 
All forms should be taken into account in any chemical model- 
ing involving deuterated ions. 

Although Pagani et al. (2009) showed that the ratio is un- 
likely in steady-state in prestellar cores, this will not affect the 
computation of the column densities since, as mentioned in 
Section 2, the collisional rates for ortho-H2D + (as well as para- 
and ortho-H2D + ) are similar for collisions with both para- and 
ortho -H2. Note also tha t the difference between our modeling 
and iPagani et al.l (|2009) comes from the time-scale used. Our 
values are at equilibrium. 

With a comparison between the chemical model predictions and 
the observations, three convergences can be found for the com- 
plete depletion case and f c / = 100. These domains are quoted in 
Table E] for densities of 10 5 , 5 x 10 5 and 1 x 10 6 cm" 3 . From 
the gas kinetic temperature range (9.5 - 12.4 K) found from the 
uncertainties on the para-D2H + and ortho-H2D + linewidths, the 
n# 2 = 10 5 cirT 3 domain can be ruled out. Should we consider the 
highest gas kinetic temperature (12.4 K) the modeling results in 
a para-D2H + /ortho-H2D + abundance ratio 12% lower than the 
observations for a molecular density of 10 5 cirT 3 . This would 
reduce the upper limit on the ortho-to-para ratio to a value of ~ 
6, ~ 40% higher than the value modeled for the complete deple- 
tion case and with a depletion factor of 100 (see Figure[5]l. Also, 
the comparison between the chemical modeling and the ortho- 
to-para D2H + upper limit from the non detection of the ortho 
transition is consistent with a low gas temperature (< 1 1.7 K for 
a 10 6 cirT 3 density). A proper model should consider the (un- 
known) physical structure and proper time-depended freeze-out. 

Note that no other molec ular tracers like DCO + and N2D + , 
observed bv lLis et aH d2002l) can be used as a comparison with 
the chemical modeling as it might prove difficult to disentangle 
their contribution from the more extended envelope to the cen- 
tral region. 

Considering an average 1 1 K gas temperature the upper limit for 
the ortho-to-para D2H + ratio is 2 to 3 times larger than the value 
found from the modeling (around 3) for densities larger than 10 5 



cm 3 . We get, from the non-LTE radiative transfer modeling, 
/V(ortho-H 2 D + ) = (1.4 + 0.1) x 10 13 cirT 2 , /V(para-D 2 H + ) = 
(1.4 ± 0.2) x 10 13 cm" 2 , and/V(ortho-D 2 H + ) < 1.4 x 10 14 cm" 2 . 
These values are reproduced by (or in the case of ortho-D2H + 
compatible with) our chemical modeling in the complete deple- 
tion case. The present state of the art instruments are clearly 
unable to detect the ortho-D2H + transition in this source. The 
CCAT (Cornell Caltech Atacama Telescop^f) project appears 
to be the most adequate to detect this ortho-D2H + transition in 
warmer regions. 

From Figure [6] we can use the average 1 1 K gas temperature in 
the complete depletion case as well as with a depletion factor 
of 100 to estimate the radiation temperature of the para-H2D + 
transition at 1370.085 GHz for densities larger than 10 5 cm -3 . 
Using a non LTE modeling described in Section 3.2, the result- 
ing para-H 2 D + column density is (9.07 ± 0.49) 10 13 cirT 2 (T R 
= (60 + 20) mK) for the complete depletion case, and (3.40 ± 
0.19) 10 13 cm" 2 (T R = (50 + 10) mK) for a depletion factor of 
100. This transition unfortunately falls in a frequency range not 
covered by the HIFI instrument, but could be a target for future 
instruments. Carbon monoxyde in the LDN 1 689N cloud is o nly 
moderately depleted in the single-dish beam dLis et al. 2002) but 
it is possible that the angular resolution of the existing single- 
dish data is simply insufficient to show the spatial stratification 
predicted by our chemical model. The presence of completely 
depleted regions smaller than the CSO beam cannot be ruled out 
by our observations. 

The deuterium peak in the LDN 1689N cloud has been argued 
in the literature to be a shock-compressed interaction region be- 
tween a molecular outflow and an ambient cloud. Therefore high 
spatial resolution mapping observations of high density tracers, 
like NH3 and its deuterated counterparts as well as N2H + and 
N2D + are necessary to allow investigating the kinematics of the 
high-density gas in this region. Both the Atacama Large Array 
Millimeter and Expanded Very Large Array are suitable for a 
follow-up of this source that reveals to be a non typical prestel- 
lar core. 



4. Conclusions 

1. The ground-state ortho-D2H + l^i-Oo.o transition at 1476.6 
GHz in the prestellar core 16293E has been searched for with 
the Herschel/HIFI instrument. 

2. The collision rates for the ortho and para H2D + and D2H + 
ions with molecular hydrogen have been used with a non- 
LTE radiative transfer code (RADEX) to derive the column 
densities of the detected ortho-H2D + and para-D2H + ground 
transitions as well as the upper limit on the ortho-D2H + ob- 
served with the HIFI instrument on board the Herschel ob- 
servatory. 

3. We used a gas-phase chemical model, in which deuterium 
chemistry and the spin states have been added and compared 
our modeling with the inferred para-D2H + /ortho-H2D + , as 
well as the upper limit on the ortho-to-para D2H + ratio. The 
kinetic temperature is deduced from the line width of both 
detected ions to be about 10 K, and we varied the molecular 
hydrogen density between 10 5 and 10 6 cm 3 . The upper limit 
on the observed ortho-to-para D2H + ratio is consistent with 
the modeling and points to a low (~ 11 K) gas kinetic tem- 
perature. The detection of the ortho-D2H + transition in the 
cold regime is a challenge as an rms of about 9 mK is needed, 



3 http://www.ccatobservatory.org/ 
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compared to the 21 mK rms reached by the HIFI instrument 
in about 5 hours. 
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Table 2. Non-LTE computations for the ortho-H2D + , para- and ortho-D2H + column densities for densities of 10 5 , 5 x 10 s and 10 6 
cm 4 . The temperature range results from the intersection of the modeling (thin lines in Figure [U and the non-LTE computation 
(thick lines in Figure @}. 



n« 2 (cm j ) 




10* 




5 x 10 s 


1 x 10" 


Tgas (K) 


[12 


.7-15.6] 


[9.9-12.2] 


[9.4-11.2] 


N(ortho-H 2 D+) 


[1.4- 


1.9] x 10" 


[1.2 


- 1.8] x 10" 


[1.2-1.9] x 10 1J 


N(para-D 2 H + ) 


[1.7- 


2.9] x 10 13 


[1.1 


- 2.4] x 10 13 


[1.1-2.3] x 10 13 


N(ortho-D 2 H + ) 


< 1 


.6x 10 14 


< 


9.2 x 10 14 


< 5.1 x 10 14 


[ortho-D 2 H + ]/[para-D 2 H+] 




< 6 




<44 


< 32 



8 



